Oscillations in turbulence-condensate system 
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We consider developed turbulence in the Gross-Pitaevsky model where condensate appears due to 
an inverse cascade. Despite being fully turbulent, the system demonstrates non-decaying periodic 
oscillations around a steady state, when turbulence and condensate periodically exchange a small 
fraction of waves. We show that these collective oscillations are not of a predator-prey type, as 
was suggested earlier; they are due to phase coherence and anomalous correlations imposed by the 
condensate. 
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Understanding the interaction of turbulence and a co- 
herent flow is an important problem in turbulence stud- 
ies in fluid mechanics and beyond, both from fundamen- 
tal and practical perspectives. In fluids, coherent flows 
are system-size vortices or zonal flows of different pro- 
flies [1, 2J, which are known to diminish the turbulence 
level, change its nature, and make its statistics more non- 
Gaussian [2 [3] . Here we consider arguably the simplest 
case of a turbulent system with a condensate: the coher- 
ent part is expected to be a constant and turbulence to 
consist of weakly interacting waves. 

The Nonlinear Schrodinger equation (NSE), also 
known as the Gross-Pitaevsky equation, provides a uni- 
versal description of the evolution of any nonlinear, spec- 
trally narrow wave packet [IHl]: 

«Vt + vV-|V'lV = o. (1) 

As such, it can provide a model of nonlinear behavior in 
a wide range of physical systems, from locally interacting 
bosons to plasmas and fluids. The addition of damping 
and driving causes the NSE to exhibit turbulence, and its 
universality makes it a major subject of interest in the 
study of wave turbulence [5-9J. Because the equation 
conserves both total energy and wave action, 

it has the potential to form both a direct and inverse 
turbulent cascade. In a system with external forcing, the 
inverse cascade leads to the growth of a single coherent 
mode known as the spectral condensate, which comes 
to dominate the dynamics of the system. Understand- 
ing the interaction between a large condensate and other 
spatial modes of a system governed by the NSE is crucial 
to a deeper understanding of wave turbulence in such a 
system. 

Numerical simulations of wave turbulence in the NSE 
suggest the existence of collective oscillations of the 
turbulence-condensate system [71 [TD] . During these oscil- 
lations, a small fraction of wave action is periodically con- 
verted from the condensate to the turbulent part of the 
spectrum, with total wave action unchanged. Such oscil- 
lations were predicted to take place when broad turbulent 



spectra coexist with a sharp spectral peak [TT] . For three- 
wave interaction, a simple model of a predator-prey type 
that describes the evolution of the total numbers of waves 
in the two groups has the form: dNo/dt — —bN^ + Nqti, 
dn I dt = — Nqu, which gives the oscillations with the 
frequency \/jb — a/ hNq [TT]. For the oscillations in 
NSE with a condensate a similar model was suggested [7] : 
dNo/dt = -bNo + Nqu'^, dn/dt = -fn - N^n. The most 
evident defect of this model is that it does not conserve 
the total number of waves. Moreover, its predictions 
for the steady-state values and the frequency of oscil- 
lations are in disagreement with the data, see [T^. It 
was also pointed out in [l^] that the account of phase 
coherence (or anomalous correlations) is needed to de- 
scribe the oscillations (for the account of anomalous cor- 
relations, see also TS"). Unfortunately, the model with 
anomalous correlations was not treated properly in |12) . 
Here we derive the simplest model of this type, based 
on anomalous correlations in the three-wave system of 
two counter-propagating waves interacting with the con- 
densate. We study the model analytically and compare 
the results with the oscillations of individual modes ex- 
tracted from the direct numerical simulations of the tur- 
bulence in the NSE. We demonstrate that the three-wave 
model provides reasonable explanation for collective os- 
cillations. 

To start, we show that an important element of turbu- 
lence against the background of a condensate must be 
anomalous correlations, i.e. phase coherence between 
waves running in opposite directions. Let us denote 
N = j-^p and Nq = The simplest condensate is 

a spatially- uniform field, ^ = ^/T7q exp{—iNQt) , which is 
an exact solution of ([T]). Small over-condensate fluctua- 
tions satisfy 

iV'fe = (fc' + 2iVo)^fe + *V-fc. (2) 

This equation gives the Bogolyubov dispersion relation, 

nl = 2Nok^ + , (3) 

for a pair of counter-propagating waves V'±fe 
exp{—iNot =F iilkt). 
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FIG. 1: (Color online) Numerical solution of ODE system (|5|-([6| with A*' = 3600 and fc = 1: phase difference (left), the 
normalized amplitude ajA — \fnfN (middle), and phase portrait (right). Different curves correspond to different initial 
conditions. 



To look into the effective dynamics of both the number 
of waves and the condensate, let us assume for simplic- 
ity that the condensate interacts with only two contra- 
propagating waves with amplitudes 

•j/j-t/j — ^/nexp^—iNot ± ikx + i4>±k) 



2nNsin9, 9 = 2k'^ + 2N{1 + cos 6*) resulting in 



and phase difference 6 = 200 



H = 2k^n + + 2n{N - 2n)(l -t- cosd) -f . (4) 

That gives the equations [Ij: 

n ^ 2n{N ~2n)suie, (5) 
9 = 2P + 2{N -An){l + cos9) + 2n. (6) 

Combining ([4| and ([5| one can obtain the equation 
for wave amplitude, fi = f{n), where f{n) = 7n^ + 



(12fc2-47V)n3^ 
AHN - 2/t2iv2^ 



{3N^ ~6H ~4:k^ ~8NP 



(AHk^ + 
The solution of 



this equation can be expressed via the elliptic functions: 
n{t) = no + If {no) {r{t, 9 2,93) - lj"{no))~\ Here 

no = + + k^- ^ is the solution of /(n) = 0, 
V is the Weierstrass elliptic function, and 92 and 93 are 
two invariants of V determined by coefficients in /(n). 

The typical solutions of system ([5])-([6]) are shown in 
Fig. [1] for a wide range of initial conditions. Motivated 
by application with very large condensate, we are mostly 
interested in limit oin/N <^ 1 (solid line in Fig. [T]). This 
limit is characterized by longer periods of oscillations, 
cusped shape of a{t) curves (where a = \/n), and by 
open trajectories in the phase space. The system spends 
most of its time around 9 — n state, avoiding both stable 
points: 6 = 0, n — {AN + fc^)/14 with its unrealistically 
high n/N ratio, and the unphysical 9 = tt, n = —\k^- 

In the limit of n ^ TV, the system ([5])-(|6| reduces to 



n{t) 
9{t) 
n{9) 



2N 



n(0) ( 1 + ^ sin^ Vtkt 



2 arctan ( tan ^Ikt 



n(0) 



27V + fc2 
A^(l + cos 61) + fc2' 



(7) 
(8) 
(9) 



Here, = V2Nk'^ + fc^ is the Bogolubov frequency, and 
n(0) = n\g^o is the constant of integration. The sec- 
ond constant of integration shifts the solution in time; it 
is selected so that 6*1 (=0 = 0. The solution ([7])-([9]) has 
the frequency of oscillations 2rifc, the stair-like time de- 
pendence of phase, and the cusped shape of amplitude, 
a = y/n oc |sinf2j.t|, as suggested by the numerical so- 
lution of shown in Fig. [l] Next, we demonstrate 
that n N approximation is well justified for the levels 
of condensate typical for NSE turbulence. 

We now compare the solution of the ODE model ([5|- 
([g]) with the results of numerical simulations of devel- 
oped turbulence in nonlinear Schrodinger equation ([T]). 
Our numerical simulations, based on the 4th order fully 




7r/3 27T/3 TT 



FIG. 2: (Color online) Spectrum of NSE turbulence nt (left) 
and phase difference 6k (right). A'' = 3600. 
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FIG. 3: (Color online) Left: oscillations in over-condensate zoomed around times t' = 800 (A'^ = 486) and t' = 3500 (A'^ = 1518). 
Middle; frequency spectrum of the time-dependence of the number of over-condensate waves for the same run in time intervals 
[800, 810] and [3500, 3510]. Right: first two frequencies as function of the number of waves next to lines with slopes 4 and 2\/2. 



dealiased split-step method [111 [15] , are set up similarly 
to [8 and are described in detail in |10) . Simulation are 
done in periodic 2tt x 27r domain. Extra terms are added 
to Eq. ([T]) to model forcing — the large scale multiplica- 
tive pumping and small scale damping. In simulations 
where thermal equilibrium is used as initial conditions, 
the pumping is needed to gradually raise the wave action 
and develop the condensate. 

In such simulations, as the wave action increases with 
time, the system undergoes the series of phase transi- 
tions (lOj. The "phases" are distinguished by different 
symmetries of the spectrum and by different spatial pat- 
terns in over-condensate fluctuations observed on small 
scales. The shape of the spectrum changes from radially 
symmetric to two-petal spectrum at ~ 500, to three- 
petal spectrum at TV ~ 1000, and to four-petal spectrum 
at ^ 2000 — 4000. Further transitions are possible. 
The transitions are related to the decrease of the angle 
of wave interactions, which causes the spectrum to break 
into a large number of narrow-angle bands. 

The typical four-petal spectrum, n^, is shown in Fig. [2] 
together with phase differences, 9k- Even though the 
phases of different modes appear to be random, the 
modes with opposite wavenumbers are pairwise corre- 
lated by their interactions with the condensate. 

Condensate oscillations at the frequencies of individ- 
ual modes add up to the overall fluctuations around con- 
densate's average, shown in Fig. [3j The highest con- 
tribution to the signal is from the modes with lowest 
wave numbers. Our three-wave model predicts that the 
modes oscillate at twice the Bogolyubov frequency. In- 
deed, the NSE simulations show that the dominant fre- 
quency of condensate-turbulence oscillations, 2^/2Nq, is 
the frequency of jkj = 1 modes, while the frequency of the 
second harmonic, 4v'iVo, correspond to jkj^ = 2 modes. 
These frequencies (which are much smaller than the fre- 
quency of the phase rotation of the condensate A^) es- 
sentially do not depend on the level of over-condensate 
fluctuations. They are clearly seen in the oscillations of 
the condensate amplitude and of the normal correlation 
functions. 



Our three-wave model describes not only the frequen- 
cies but also the time dependence of the oscillations. The 
time dependence of the phase and the amplitude of an 
individual mode extracted from a NSE simulation is com- 
pared to the prediction of the model ([5|-([6| in Fig. B] 
Here, the n/N ratio is small, and the reduced model ^j- 
^ works as good as the full one, as illustrated in phase 
portrait section of Fig. |4] As predicted by the model, 
the amplitude has cusped shape, while the phase is 
localized around 6 = tt. 

The biggest difference between the NSE modes and 
the model is the phase portrait. The model predicts 
monotonous increase of the phase from (2j — l)7r to 
(2j + l)7r, while the phase of the NSE modes oscillates 
around ±7r in closed loops. Apparently, the interac- 
tion with other modes, unaccounted in our simple model, 
leads to phase locking. 

The amplitude in the three-wave model is determined 
up to the constant of integration. In other word, the 
level of fluctuations must be obtained from the turbu- 
lence data. Our study of NSE turbulence at different 
levels of condensate indicates that the amplitude of lower 
modes remain roughly constant for the wide range of A'^. 
This also means that n/N decreases as the condensate 
level increases, as shown in Fig. [5] (There is a possibility 
that n/N is somewhat larger near phase transition but 
this topic needs additional investigation.) In all cases 
considered, the ratio n/N is small, well within the appli- 
cability limits of the reduced model. At higher A^, the 
loops in the phase portrait are reconnecting tighter and 
closer to 6* = TT. It is also interesting that amplitudes of 
non-condensate modes are larger than fluctuations of the 
condensate. 

Regarding the data presented in Fig. [5] we shall make 
the following technical comment, unaddressed in [lOj . 
The states, similar to the state shown in Fig. [2] can be 
achieved in numerical simulation by using initial condi- 
tions where the condensate is superimposed on top of 
the thermal equilibrium spectrum. For such systems, 
it take only a few linear time units to relax to quasi- 
equilibrium with the appropriate symmetry, in contrast 
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FIG. 4: (Color online) First three panels: comparison of oscillations of k = (1, 1) mode extracted from simulation of turbulence 
at A*' = 400 with the prediction of three-wave model; from left to right, phase, amplitude, the phase portrait are shown. Fourth 
panel: comparison of phase portraits for A'' = 3600 and k = (0,3). Here, as in Fig. 1 ajA = ^n/N . 
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FIG. 5: (Color online) Amplitude of oscillations of low modes 
for different levels of condensate is averaged over time interval 
[10, 12] in simulations with preset condensate. 



to thousands of time units of evolution needed for sys- 
tems without preset condensate. Moreover, as shown in 



Fig. [3| the condensate-turbulence oscillations are very 
fast and can be studied in simulations without forcing 
(except for the small damping for smooth transition to 
de-aliased region). 

To summarize, at the large level of condensate the 
three-wave model capture the following features of the 
turbulence-condensate oscillations: (i) the frequency of 
oscillations is twice the Bogolubov frequency, (ii) the sys- 
tem spends most of its time around 9 — n state, and (iii) 
the amplitude as a function of time has a non-trivial, 
cusped shape. To describe the shape of the phase os- 
cillations around ±7r, additional mechanisms need to be 
included. 
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